The neuroendocrine transition in prostate cancer is dynamic and dependent on ASCL1.

Lineage plasticity is a recognized hallmark of cancer progression that can shape therapy outcomes. The underlying cellular and molecular mechanisms mediating lineage plasticity remain poorly understood. Here, we describe a versatile in vivo platform to identify and interrogate the molecular determinants of neuroendocrine lineage transformation at different stages of prostate cancer progression. Adenocarcinomas reliably develop following orthotopic transplantation of primary mouse prostate organoids acutely engineered with human-relevant driver alterations (e.g., Rb1-/-; Trp53-/-; cMyc+ or Pten-/-; Trp53-/-; cMyc+), but only those with Rb1 deletion progress to ASCL1+ neuroendocrine prostate cancer (NEPC), a highly aggressive, androgen receptor signaling inhibitor (ARSI)-resistant tumor. Importantly, we show this lineage transition requires a native in vivo microenvironment not replicated by conventional organoid culture. By integrating multiplexed immunofluorescence, spatial transcriptomics and PrismSpot to identify cell type-specific spatial gene modules, we reveal that ASCL1+ cells arise from KRT8+ luminal epithelial cells that progressively acquire transcriptional heterogeneity, producing large ASCL1+;KRT8- NEPC clusters. Ascl1 loss in established NEPC results in transient tumor regression followed by recurrence; however, Ascl1 deletion prior to transplantation completely abrogates lineage plasticity, yielding adenocarcinomas with elevated AR expression and marked sensitivity to castration. The dynamic feature of this model reveals the importance of timing of therapies focused on lineage plasticity and offers a platform for identification of additional lineage plasticity drivers.

Metastatic site influences driver gene function in pancreatic cancer.

Driver gene mutations can increase the metastatic potential of the primary tumor1-3, but their role in sustaining tumor growth at metastatic sites is poorly understood. A paradigm of such mutations is inactivation of SMAD4 - a transcriptional effector of TGFß signaling - which is a hallmark of multiple gastrointestinal malignancies4,5. SMAD4 inactivation mediates TGFß's remarkable anti- to pro-tumorigenic switch during cancer progression and can thus influence both tumor initiation and metastasis6-14. To determine whether metastatic tumors remain dependent on SMAD4 inactivation, we developed a mouse model of pancreatic ductal adenocarcinoma (PDAC) that enables Smad4 depletion in the pre-malignant pancreas and subsequent Smad4 reactivation in established metastases. As expected, Smad4 inactivation facilitated the formation of primary tumors that eventually colonized the liver and lungs. By contrast, Smad4 reactivation in metastatic disease had strikingly opposite effects depending on the tumor's organ of residence: suppression of liver metastases and promotion of lung metastases. Integrative multiomic analysis revealed organ-specific differences in the tumor cells' epigenomic state, whereby the liver and lungs harbored chromatin programs respectively dominated by the KLF and RUNX developmental transcription factors, with Klf4 depletion being sufficient to reverse Smad4's tumor-suppressive activity in liver metastases. Our results show how epigenetic states favored by the organ of residence can influence the function of driver genes in metastatic tumors. This organ-specific gene-chromatin interplay invites consideration of anatomical site in the interpretation of tumor genetics, with implications for the therapeutic targeting of metastatic disease.

Microbial cancer immunotherapy reprograms hematopoietic stem cells to enhance anti-tumor immunity.

Mycobacterium bovis BCG is the vaccine against tuberculosis and an immunotherapy for bladder cancer. When administered intravenously, BCG reprograms bone marrow hematopoietic stem and progenitor cells (HSPCs), leading to heterologous protection against infections. Whether HSPC-reprogramming contributes to the anti-tumor effects of BCG administered into the bladder is unknown. We demonstrate that BCG administered in the bladder in both mice and humans reprograms HSPCs to amplify myelopoiesis and functionally enhance myeloid cell antigen presentation pathways. Reconstitution of naive mice with HSPCs from bladder BCG-treated mice enhances anti-tumor immunity and tumor control, increases intratumoral dendritic cell infiltration, and synergizes with checkpoint blockade. We conclude that bladder BCG acts systemically, reprogramming HSPC-encoded innate immunity, highlighting the broad potential of modulating HSPC phenotypes to improve tumor immunity. One Sentence Summary: BCG administered in the bladder reprograms bone marrow HSPCs and contributes to tumor control via enhanced myeloid cells.

The covariance environment defines cellular niches for spatial inference.

A key challenge of analyzing data from high-resolution spatial profiling technologies is to suitably represent the features of cellular neighborhoods or niches. Here we introduce the covariance environment (COVET), a representation that leverages the gene-gene covariate structure across cells in the niche to capture the multivariate nature of cellular interactions within it. We define a principled optimal transport-based distance metric between COVET niches that scales to millions of cells. Using COVET to encode spatial context, we developed environmental variational inference (ENVI), a conditional variational autoencoder that jointly embeds spatial and single-cell RNA sequencing data into a latent space. ENVI includes two decoders: one to impute gene expression across the spatial modality and a second to project spatial information onto single-cell data. ENVI can confer spatial context to genomics data from single dissociated cells and outperforms alternatives for imputing gene expression on diverse spatial datasets.

Novel Approach for Pancreas Transcriptomics Reveals the Cellular Landscape in Homeostasis and Acute Pancreatitis.

BACKGROUND & AIMS: Acinar cells produce digestive enzymes that impede transcriptomic characterization of the exocrine pancreas. Thus, single-cell RNA-sequencing studies of the pancreas underrepresent acinar cells relative to histological expectations, and a robust approach to capture pancreatic cell responses in disease states is needed. We sought to innovate a method that overcomes these challenges to accelerate study of the pancreas in health and disease. METHODS: We leverage FixNCut, a single-cell RNA-sequencing approach in which tissue is reversibly fixed with dithiobis(succinimidyl propionate) before dissociation and single-cell preparation. We apply FixNCut to an established mouse model of acute pancreatitis, validate findings using GeoMx whole transcriptome atlas profiling, and integrate our data with prior studies to compare our method in both mouse and human pancreas datasets. RESULTS: FixNCut achieves unprecedented definition of challenging pancreatic cells, including acinar and immune populations in homeostasis and acute pancreatitis, and identifies changes in all major cell types during injury and recovery. We define the acinar transcriptome during homeostasis and acinar-to-ductal metaplasia and establish a unique gene set to measure deviation from normal acinar identity. We characterize pancreatic immune cells, and analysis of T-cell subsets reveals a polarization of the homeostatic pancreas toward type-2 immunity. We report immune responses during acute pancreatitis and recovery, including early neutrophil infiltration, expansion of dendritic cell subsets, and a substantial shift in the transcriptome of macrophages due to both resident macrophage activation and monocyte infiltration. CONCLUSIONS: FixNCut preserves pancreatic transcriptomes to uncover novel cell states during homeostasis and following pancreatitis, establishing a broadly applicable approach and reference atlas for study of pancreas biology and disease.

Deep generative model deciphers derailed trajectories in acute myeloid leukemia.

Single-cell genomics has the potential to map cell states and their dynamics in an unbiased way in response to perturbations like disease. However, elucidating the cell-state transitions from healthy to disease requires analyzing data from perturbed samples jointly with unperturbed reference samples. Existing methods for integrating and jointly visualizing single-cell datasets from distinct contexts tend to remove key biological differences or do not correctly harmonize shared mechanisms. We present Decipher, a model that combines variational autoencoders with deep exponential families to reconstruct derailed trajectories (https://github.com/azizilab/decipher). Decipher jointly represents normal and perturbed single-cell RNA-seq datasets, revealing shared and disrupted dynamics. It further introduces a novel approach to visualize data, without the need for methods such as UMAP or TSNE. We demonstrate Decipher on data from acute myeloid leukemia patient bone marrow specimens, showing that it successfully characterizes the divergence from normal hematopoiesis and identifies transcriptional programs that become disrupted in each patient when they acquire NPM1 driver mutations.

Supervised discovery of interpretable gene programs from single-cell data.

Factor analysis decomposes single-cell gene expression data into a minimal set of gene programs that correspond to processes executed by cells in a sample. However, matrix factorization methods are prone to technical artifacts and poor factor interpretability. We address these concerns with Spectra, an algorithm that combines user-provided gene programs with the detection of novel programs that together best explain expression covariation. Spectra incorporates existing gene sets and cell-type labels as prior biological information, explicitly models cell type and represents input gene sets as a gene-gene knowledge graph using a penalty function to guide factorization toward the input graph. We show that Spectra outperforms existing approaches in challenging tumor immune contexts, as it finds factors that change under immune checkpoint therapy, disentangles the highly correlated features of CD8+ T cell tumor reactivity and exhaustion, finds a program that explains continuous macrophage state changes under therapy and identifies cell-type-specific immune metabolic programs.

scKINETICS: inference of regulatory velocity with single-cell transcriptomics data.

MOTIVATION: Transcriptional dynamics are governed by the action of regulatory proteins and are fundamental to systems ranging from normal development to disease. RNA velocity methods for tracking phenotypic dynamics ignore information on the regulatory drivers of gene expression variability through time. RESULTS: We introduce scKINETICS (Key regulatory Interaction NETwork for Inferring Cell Speed), a dynamical model of gene expression change which is fit with the simultaneous learning of per-cell transcriptional velocities and a governing gene regulatory network. Fitting is accomplished through an expectation-maximization approach designed to learn the impact of each regulator on its target genes, leveraging biologically motivated priors from epigenetic data, gene-gene coexpression, and constraints on cells' future states imposed by the phenotypic manifold. Applying this approach to an acute pancreatitis dataset recapitulates a well-studied axis of acinar-to-ductal transdifferentiation whilst proposing novel regulators of this process, including factors with previously appreciated roles in driving pancreatic tumorigenesis. In benchmarking experiments, we show that scKINETICS successfully extends and improves existing velocity approaches to generate interpretable, mechanistic models of gene regulatory dynamics. AVAILABILITY AND IMPLEMENTATION: All python code and an accompanying Jupyter notebook with demonstrations are available at http://github.com/dpeerlab/scKINETICS.

The covariance environment defines cellular niches for spatial inference.

The tsunami of new multiplexed spatial profiling technologies has opened a range of computational challenges focused on leveraging these powerful data for biological discovery. A key challenge underlying computation is a suitable representation for features of cellular niches. Here, we develop the covariance environment (COVET), a representation that can capture the rich, continuous multivariate nature of cellular niches by capturing the gene-gene covariate structure across cells in the niche, which can reflect the cell-cell communication between them. We define a principled optimal transport-based distance metric between COVET niches and develop a computationally efficient approximation to this metric that can scale to millions of cells. Using COVET to encode spatial context, we develop environmental variational inference (ENVI), a conditional variational autoencoder that jointly embeds spatial and single-cell RNA-seq data into a latent space. Two distinct decoders either impute gene expression across spatial modality, or project spatial information onto dissociated single-cell data. We show that ENVI is not only superior in the imputation of gene expression but is also able to infer spatial context to disassociated single-cell genomics data.

Conserved transcriptional connectivity of regulatory T cells in the tumor microenvironment informs new combination cancer therapy strategies.

While regulatory T (Treg) cells are traditionally viewed as professional suppressors of antigen presenting cells and effector T cells in both autoimmunity and cancer, recent findings of distinct Treg cell functions in tissue maintenance suggest that their regulatory purview extends to a wider range of cells and is broader than previously assumed. To elucidate tumoral Treg cell 'connectivity' to diverse tumor-supporting accessory cell types, we explored immediate early changes in their single-cell transcriptomes upon punctual Treg cell depletion in experimental lung cancer and injury-induced inflammation. Before any notable T cell activation and inflammation, fibroblasts, endothelial and myeloid cells exhibited pronounced changes in their gene expression in both cancer and injury settings. Factor analysis revealed shared Treg cell-dependent gene programs, foremost, prominent upregulation of VEGF and CCR2 signaling-related genes upon Treg cell deprivation in either setting, as well as in Treg cell-poor versus Treg cell-rich human lung adenocarcinomas. Accordingly, punctual Treg cell depletion combined with short-term VEGF blockade showed markedly improved control of PD-1 blockade-resistant lung adenocarcinoma progression in mice compared to the corresponding monotherapies, highlighting a promising factor-based querying approach to elucidating new rational combination treatments of solid organ cancers.

Epigenetic plasticity cooperates with cell-cell interactions to direct pancreatic tumorigenesis.

The response to tumor-initiating inflammatory and genetic insults can vary among morphologically indistinguishable cells, suggesting as yet uncharacterized roles for epigenetic plasticity during early neoplasia. To investigate the origins and impact of such plasticity, we performed single-cell analyses on normal, inflamed, premalignant, and malignant tissues in autochthonous models of pancreatic cancer. We reproducibly identified heterogeneous cell states that are primed for diverse, late-emerging neoplastic fates and linked these to chromatin remodeling at cell-cell communication loci. Using an inference approach, we revealed signaling gene modules and tissue-level cross-talk, including a neoplasia-driving feedback loop between discrete epithelial and immune cell populations that was functionally validated in mice. Our results uncover a neoplasia-specific tissue-remodeling program that may be exploited for pancreatic cancer interception.

Cross-platform dataset of multiplex fluorescent cellular object image annotations.

Defining cellular and subcellular structures in images, referred to as cell segmentation, is an outstanding obstacle to scalable single-cell analysis of multiplex imaging data. While advances in machine learning-based segmentation have led to potentially robust solutions, such algorithms typically rely on large amounts of example annotations, known as training data. Datasets consisting of annotations which are thoroughly assessed for quality are rarely released to the public. As a result, there is a lack of widely available, annotated data suitable for benchmarking and algorithm development. To address this unmet need, we release 105,774 primarily oncological cellular annotations concentrating on tumor and immune cells using over 40 antibody markers spanning three fluorescent imaging platforms, over a dozen tissue types and across various cellular morphologies. We use readily available annotation techniques to provide a modifiable community data set with the goal of advancing cellular segmentation for the greater imaging community.

Modeling collective cell behavior in cancer: Perspectives from an interdisciplinary conversation.

Collective cell behavior contributes to all stages of cancer progression. Understanding how collective behavior emerges through cell-cell interactions and decision-making will advance our understanding of cancer biology and provide new therapeutic approaches. Here, we summarize an interdisciplinary discussion on multicellular behavior in cancer, draw lessons from other scientific disciplines, and identify future directions.

Distinct tumor architectures for metastatic colonization of the brain.

Brain metastasis is a dismal cancer complication, hinging on the initial survival and outgrowth of disseminated cancer cells. To understand these crucial early stages of colonization, we investigated two prevalent sources of cerebral relapse, triple-negative (TNBC) and HER2+ breast cancer (HER2BC). We show that these tumor types colonize the brain aggressively, yet with distinct tumor architectures, stromal interfaces, and autocrine growth programs. TNBC forms perivascular sheaths with diffusive contact with astrocytes and microglia. In contrast, HER2BC forms compact spheroids prompted by autonomous extracellular matrix components and segregating stromal cells to their periphery. Single-cell transcriptomic dissection reveals canonical Alzheimer's disease-associated microglia (DAM) responses. Differential engagement of tumor-DAM signaling through the receptor AXL suggests specific pro-metastatic functions of the tumor architecture in both TNBC perivascular and HER2BC spheroidal colonies. The distinct spatial features of these two highly efficient modes of brain colonization have relevance for leveraging the stroma to treat brain metastasis.

Single-cell analysis of bidirectional reprogramming between early embryonic states reveals mechanisms of differential lineage plasticities.

Two distinct fates, pluripotent epiblast (EPI) and primitive (extra-embryonic) endoderm (PrE), arise from common progenitor cells, the inner cell mass (ICM), in mammalian embryos. To study how these sister identities are forged, we leveraged embryonic (ES) and eXtraembryonic ENdoderm (XEN) stem cells - in vitro counterparts of the EPI and PrE. Bidirectional reprogramming between ES and XEN coupled with single-cell RNA and ATAC-seq analyses uncovered distinct rates, efficiencies and trajectories of state conversions, identifying drivers and roadblocks of reciprocal conversions. While GATA4-mediated ES-to-iXEN conversion was rapid and nearly deterministic, OCT4, KLF4 and SOX2-induced XEN-to-iPS reprogramming progressed with diminished efficiency and kinetics. The dominant PrE transcriptional program, safeguarded by Gata4, and globally elevated chromatin accessibility of EPI underscored the differential plasticities of the two states. Mapping in vitro trajectories to embryos revealed reprogramming in either direction tracked along, and toggled between, EPI and PrE in vivo states without transitioning through the ICM.

SEACells infers transcriptional and epigenomic cellular states from single-cell genomics data.

Metacells are cell groupings derived from single-cell sequencing data that represent highly granular, distinct cell states. Here we present single-cell aggregation of cell states (SEACells), an algorithm for identifying metacells that overcome the sparsity of single-cell data while retaining heterogeneity obscured by traditional cell clustering. SEACells outperforms existing algorithms in identifying comprehensive, compact and well-separated metacells in both RNA and assay for transposase-accessible chromatin (ATAC) modalities across datasets with discrete cell types and continuous trajectories. We demonstrate the use of SEACells to improve gene-peak associations, compute ATAC gene scores and infer the activities of critical regulators during differentiation. Metacell-level analysis scales to large datasets and is particularly well suited for patient cohorts, where per-patient aggregation provides more robust units for data integration. We use our metacells to reveal expression dynamics and gradual reconfiguration of the chromatin landscape during hematopoietic differentiation and to uniquely identify CD4 T cell differentiation and activation states associated with disease onset and severity in a Coronavirus Disease 2019 (COVID-19) patient cohort.

STING inhibits the reactivation of dormant metastasis in lung adenocarcinoma.

Metastasis frequently develops from disseminated cancer cells that remain dormant after the apparently successful treatment of a primary tumour. These cells fluctuate between an immune-evasive quiescent state and a proliferative state liable to immune-mediated elimination1-6. Little is known about the clearing of reawakened metastatic cells and how this process could be therapeutically activated to eliminate residual disease in patients. Here we use models of indolent lung adenocarcinoma metastasis to identify cancer cell-intrinsic determinants of immune reactivity during exit from dormancy. Genetic screens of tumour-intrinsic immune regulators identified the stimulator of interferon genes (STING) pathway as a suppressor of metastatic outbreak. STING activity increases in metastatic progenitors that re-enter the cell cycle and is dampened by hypermethylation of the STING promoter and enhancer in breakthrough metastases or by chromatin repression in cells re-entering dormancy in response to TGFß. STING expression in cancer cells derived from spontaneous metastases suppresses their outgrowth. Systemic treatment of mice with STING agonists eliminates dormant metastasis and prevents spontaneous outbreaks in a T cell- and natural killer cell-dependent manner-these effects require cancer cell STING function. Thus, STING provides a checkpoint against the progression of dormant metastasis and a therapeutically actionable strategy for the prevention of disease relapse.

Leptomeningeal anti-tumor immunity follows unique signaling principles.

Metastasis to the cerebrospinal fluid (CSF)-filled leptomeninges, or leptomeningeal metastasis (LM), represents a fatal complication of cancer. Proteomic and transcriptomic analyses of human CSF reveal a substantial inflammatory infiltrate in LM. We find the solute and immune composition of CSF in the setting of LM changes dramatically, with notable enrichment in IFN-γ signaling. To investigate the mechanistic relationships between immune cell signaling and cancer cells within the leptomeninges, we developed syngeneic lung, breast, and melanoma LM mouse models. Here we show that transgenic host mice, lacking IFN-γ or its receptor, fail to control LM growth. Overexpression of Ifng through a targeted AAV system controls cancer cell growth independent of adaptive immunity. Instead, leptomeningeal IFN-γ actively recruits and activates peripheral myeloid cells, generating a diverse spectrum of dendritic cell subsets. These migratory, CCR7+ dendritic cells orchestrate the influx, proliferation, and cytotoxic action of natural killer cells to control cancer cell growth in the leptomeninges. This work uncovers leptomeningeal-specific IFN-γ signaling and suggests a novel immune-therapeutic approach against tumors within this space.